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It is well known that Hermite polynomial (HP) is an important special function, which forms a class of orthogonal 
and complete set in function space[l|-Q. Among them, the two- variable Hermite polynomial (TVHP) -ff m , n (0 0)[IIIII 
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In this paper, by virtue of the entangled state representation we concisely derive some new oper- 
ator identities regarding to two- variable Hermite polynomial (TVHP). By them and the technique 
of integration within an ordered product (IWOP) of operators we further derive new generating 
function formulas of TVHP. They are useful in quantum optical theoretical calculations. It is seen 
from this work that by combining the IWOP technique and quantum mechanical representations 
one can derive some new integration formulas even without really performing the integration. 
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> ! I. INTRODUCTION 

o 

, ; mi ^ >n] (-iy m 'ri'e™-<e n ^ 

flm,»(£,0= L H(m-l)\(n-l)\ (1) 

has become more and more useful (note that -ff mj n(0 O) is not a direct product of two independent single- variable HP). 
For examples, in quantum optics theory if m , n (0 O) can be considered as the basis of generalized Bargmann function 
space corresponding to the two-mode Fock state[5(, and is very useful for discussing quasi-probability distribution 
of some quantum states[6[. In Fourier optics theory, -Hm,n(0O) is proved to be the eigenfunction of the complex 
fractional Fourier transform, the corresponding phenomenon may be observed in the light propagation in graded index 
, (GRIN) medium[7|, this eigenmode also exists in two-dimensional Talbot effect demonstrated in GRIN medium|8]. 
J> ■ In Ref.0, it is pointed out that the squeezed two-mode number state is just a TVHP excitation on the two-mode 
[ squeezed vacuum state. 

Recalling that the entangled state representation |£) is expressed asfTTlj 

\Z) = e -m a +&W-*W\m), £ = ei+if2) (2) 

which obeys 

(a + ^10=00, (« f +b) |0= 010, (3) 



where 1 00) is the two-mode vacuum state, a' and are the Bose creation operator with [a,a^\ = [b,b<\ = 1. |0 is 
orthonormal and complete, 

fl |f) (0 = [ H, e -{a"+b-e)(a+^-i) . = 1; (4) 

where : : stands for normal ordering, and we have used the technique of integration within an ordered product 
(IWOP) of operators [l2l - tl6l | as well as |00) (00 1 =: exp[— a^a — b^b] : . According to the generating function of 
H m ,n (00) 



E 



fin f/n 

t -^H mtn (t,e) = e- tt ' + * +t 't\ (5) 
— m\n\ 

m,n— 

we can expand 

I0 = e-^'«i 2 ^—rr H m,n (OO)loo), (6) 

z — ' mini 

m,n—0 
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then the overlap between (£| and the two-mode Fock state is 

<£| m, n) = e -^=H^ n (& n , |m, n> = |00> . (7) 

Eq.([6|) exhibits quantum entanglement (l7| . In this work we shall employ the entangled state representation for deriving 
some new operator identities regarding to H m ^ n (£, £*) and then present their applications. By virtue of them, we can 
also derive some new generating function formulas of H m .n (£, £*) , which is quite useful in calculating normalization of 
some quantum states. We conclude that by combining the IWOP technique and quantum mechanical representations 
one can derive some new integration formulas even without really performing the integration. 

II. NEW OPERATOR IDENTITIES REGARDING TO TWO- VARIABLE HERMITE POLYNOMIALS 

In order to obtain the normally ordered expansion of H m n (a + b\ + fe) , by noticing [a + , + b\ = 0, we 
have 

-tt'+t(a+b r ')+t'(a 1 <+b') t(a+6 t )+i'(a t +&) . 



00 j-m^/n 



mini 

m,n— 



where we have used Baker-Hausdorff formula 

e A e B = e B e A e [A ' B] , [A, [A, B]] = [B, [A, B]] = 0. (9) 
Comparing with Eq.® we see the identity 

fl TO ,„(o + 6 t ,o t + 6) =: (a + 6 t ) m (a t +6)" : . (10) 
Using Eqs.© and © as well as the IWOP technique [I3 - [l§ |. we have 



H m , n {a + b\a) +b) = I ^({.OIO^I 



7T 



= /^ff m) n«,r):e-( at + 6 -^)(«+ 6t -0: 

= : (a + 6t) m ( a t +6 )" : , (ii) 



which implies an integration formula 

-;2 



J ^H m , n (£, D exp [- (a* - O (« " 0] = a™"*™- ( 12 ) 

Thus we see that by combining the IWOP technique and quantum mechanical representations one can derive some 
new integration formulas even without really performing the integration. Moreover, from the antinormally ordered 

expansion (denoted by : :) 



e -tt'+t'(a t +f))+t(a+6 t ) _ : e -2tt'+t'(a t +b)+t(a+6 t ): 

v o(m+n) /2-i-rti-i-ln / n 

\T - 1 -u _ [ _ 

V2 ' V2 



°^ 2 (m+n)/2 t m t >n . / a + &t a t + &\ . 

2- — — '- Hmn [—FT'—fr )'■ (13) 



m,n— 

and comparing Eq. ((5J with Eq. ()13j) we see 



H m , n (a + fct, t + 6 ) = 2t-+«)/ 2 ; j ff m>n , ^) 



(14) 
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On the other hand, it is seen that 



00 ^m-^ln 



mm 

m,n— 



_ e t(a+tf) e t' {at +b) 

__ t*'+t'( fc + at ) e t ( a + bt ) • 

= ; e - i * i *' + t '( f '+ at ) e f ( a + f ' t ) ; (15) 

= V 1 j , l , J : ff m ,„ [-i (a + fot) , _< ( a t + &)] : , (16) 
* — ' mini 

m,n— 

this leads to the new operator identity 

(a + &t) m ( a t + b ) n = im+n . Hm n [_ i ( fl + 6 t) i _ 4 ( a t + 6)] : . (17) 

Similarly, using Eqs.Q and ([4]), we also have 

(a+6t) m ( a t + & )" = / ^rr n ie) (ci 

= / r" : exp [- (at + 6 _ £*) ( a + 6 t _ f )] . 

= i m+n : tf m ,„ [-i (a + fet) , ( a t +b)]:, (18) 

which implies another integration formula 

J ^A^^n exp [_ {a * _ e) {a _ c)] = ^+n Hm n ( _ iaj _ ia *) _ (19) 

This is the reciprocal relation of Eq. (fT2j) . It is clear shown that Eqs. (fT2]) and (fl~9| are the mutual integration 
transformation between i? m ,„ (£,£*) an d the polynomial £ m £*™. 

To derive another more complicated generating function formula about TVHP, we introduce another two-mode 
entangled state involved in modes c\d\ 

10 = e-^l'+^+f*' 11 -^ |00) , C = Ci + i£>. (20) 
where \c, ct] = 1 and [d, tfn = 1, obeying 

(c + rf t )|C)=CIC), ( C t + d)|C)=C*IC), (21) 

and 



^| C )(C|= /^ :e -(^-C)(c + ^-C) ; = i. (22) 

7T ./ 7T 



As a result of Eq. (|TT|) we deduce 



— gin 

G = V - 1 —rH m , n (a + &t, at + 6) iT m ,„ (c + tft c t + rf) 

m,n— 
m,n— 



In reference to 



exp [s (a + &t) ( c + dt) +t(a* + b) (c f + d)] : . (23) 



i± eV \z\*+fz+gz* = _i e -fg/v : Rer? <0, (24) 
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the right-hand side of Eq. (1231 can be the result of the following integration with use of the IWOP technique 



. e -|C| 3 +C[(c + +d)+-(a+6+)]+C'[(H-d t )+t(« t +6)] 
x e -st(a i +b)(a+b i )~(c i +d)(c+d i ) . 

— 10 (CI e~ st ( at+b )( a+bt ) +sC ( a+,,t ) +t ' ; *( at+^, ) • 



G 



(25) 

where, within the normal ordering : :, the terms involving (a + i/M and (a) + 6) can be expressed as the following 
integration 



-st(a t +&)(a+b t )+sC(a+& t )+tC*(a t +f)) 



1 IT 6XP 









1 - is 



+r ( a + 6 t + _^_( a t + 6 ) ( a + 6 t) 



1 - ts) 1 - ts 

Substituting Eq.fl2B]) into Ea.([2"5]) and using Eqs.© and (gj| as well as Eqs.fl2TJ and we see 



(26) 



G 



1 fd 2 e fd 2 ( 



x exp 



, , , , 10 ® 10 (CI ® (CI 

^(ici 2 + ici 2 ) , <c , *c*e* 



1 - ts 



1 - ts 1 - ts 



1 



■ exp 



-ts 



— [( a t + 6) ( a + 6t) + ( c t + rf) ( c + d t)] 



+ 



1 - ts 

s (a + 6+) (c + d+) t ( c t + d) ( a t + b) 



1 - ts 



1 - ts 



(27) 



Since (a* + 6) , (a + irj j ( c t _|_ anc i ( c gjtj are a n commutative among themselves, we can make replacement 
(a 1 * + b) -> y, (a + b t ) -> x, (c* + d) -> y', (c + d*) -> x', then by comparing Ea.([2"7P with Eq.([23J we obtain the 
complicated generating function formula about TVHP 



— — H m>n (x, y) H m n (x', y') = — exp 



m\n 



1 - ts 



sxa/ + tyy' — ts (xy + x'y') 



1 - ts 



(28) 



III. APPLICATIONS 



We now present some applications of Eq. (|28|) . To begin with, we point out that Eq. ([28|) can be used for deriving 
another generating function formula. In fact, using the generating function of Laguerre polynomial^ 



^L m ( a; ) s - = (l- s )- 1 e X p( T f^j, (29) 

m=0 ^ ' 

and its relation to TVHP 

H mi m (x, y) = (-1)™ m\L m {xy) , (30) 



we can reexpress Eq. f|28|) as 



girij/n, 

(x',y') 

£ — ' m\n\ 

m^n—O 



e tyy 
1-st 



exp 



— (xy + x'y' — tyy' — \xx'\ ts 



1 - st 



771—0 

OO 



e^' ( st ) m L ™ Uy + x'y 1 - tyy' - - 



= e tvv' J2 { -^-H n 
m=0 



rr> 



i(Viy' - —=), i{\fty - —= 



Comparing the same power of s on the above two sides yields 

00 t n r x x' 

^ —\ H m,n (x,y)H„ hn (x',y) = (-t) m e tyy H m>m i(Vt,y' - -j=),i(Viy y=) 



n=0 

Further, by noticing 



gt'agiat _ gtt'gta+gt'a _. g ( — it') ia+( — it )ia 1 " — (— it) ( — it' ) . 

= > r-j : -ff m ,n (ia T ,ia) : 

' m!n! v ' 

m,n— 



and comparing it with 



e f'a e ta+ 



E 



n!m! 



a"a tm , 



we have the compact operator identity 



i n a^ = {-i) m+n :H m<n (\a\ia) 



It then follows from Ed.([28]l that 



„sa& to" b^ 

e e 



m=Q n— 

°° (— \y n+n s m t n 



m,n—0 

l 

1 - ts 



exp 



(at a + 6+6) + sa t 6 t + ia6 



1 - ts 



_ 1 c -rj TT a t b t c -fa t a+b t b) ln(l-ts) gi^ab 

1 - is 

where at the last step we have used[l5j 

e Aata =: exp[(e A -l)a t a]: . 

This result of Eq. (|5rj|) agree with that of Ref.[? ]. On the other hand, using the antinormally ordered 
■H mt „ (a\ a) : = a^ m a n , where using its P-representation as well as Eq. (TT2)) we easily prove that 



H m . n (a\a) i = / — H m>n (z*,z) \z) (z\ = I —H m>n (z*,z) : e-^-**)*-*) : = J m a n 



7T 
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we can obtain 



J-TL TTL 

J2 ^-ra t "a m 6 t "6 m 



7 1 m I 



m,n— 



E ^ 



m,n— 



1 - is 



m!n! 



(a,a t ) /f m ,„ (6, & f ) : 



exp 



-ts (a^a + tfb) + ta^tf + sab 



1 - ts 



1 - ts 



(39) 



Finally, according to Eqs.(f30|). (|32|) and (f3"5j) . we also have 



mjm a* 6^ tanhA 

a e 



E tanh t „ ^ tw 

71.1 



ri=0 



(-ir E ( tan , hA) " : ^» ( ia - iat ) *w ( i6 < i6t ) 



n=0 

m!tanh m Ae otfet tanhA : 



L m (— a^a — 6^6 — a& coth A — a^V tanh A) 



(40) 



which relates to Laguerre polynomial. Based on this, we obtain that the two-mode photon-subtracted squeezed 
vacuum state is expressed as[19f 



|A) m = a m 6 m S 2 (A)|00) 

m 7 m tanh A 

= a o e 



|00) 



= to! tanh" 1 Ae atbt tanhA L m (-a t 6 t tanhA) |00) 
= ml tanh" 1 AL m (-aW tanhA) 5 2 (A) |00) , 



(41) 



where ^(A) = exp [A(a^6^ — a&)l is the two-mode squeezing operator with A being a real squeezing parameter. From 
Eq. (|4ip . | A) can is equivalent to Laguerre polynomial excitation on squeezed vacuum state. Recall in Ref . [20| . we 
have calculated its normalization factor 

(00| e ^tanh V m & tm a m 6 m e atf,t tanhA |00) = (m!) 2 sinh 2m XP m (cosh2A) (42) 

where P m (x) is Legendre polynomial [J [l9| 

[to/2] i ^ _ 2^)!a; m " 2/ 

Pm(:E) = £ 2"/!(m-0!(m-2Z)!- (43) 

Using Eqs. (|4*Tj) and (jUJ) as well as the coherent state's completeness relation J d 13 \a,0) (a, /3\ = 1, it follows 

d2a f P L m (-a/3 tanh A) L m (-a* (3* tanh A) e -M 2 -l/3| 2 +M+a*r ) tanh a = cosh 2,„ ^ ( cogh 2A j ^ 



/ 



7T- 



which is a new integration formula. 

In summary, by virtue of the entangled state representation we concisely derive some new operator identities 
regarding to two-variable Hermite polynomials. They are useful in quantum optical theoretical calculations. By 
combining the IWOP technique and quantum mechanical representations one can derive some new integration formulas 
even without really performing the integration. 
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